Generalized fluctuation relation and effective temperatures in a driven fluid 
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By numerical simulation of a Lennard- Jones like liquid driven by a velocity gradient 7 we test 
the fluctuation relation (FR) below the (numerical) glass transition temperature Tg. We show 
that, in this region, the FR deserves to be generalized introducing a numerical factor X(r, 7) < 1 
that defines an "effective temperature" Tfr = T/X. On the same system we also measure the 
effective temperature T^ff, as defined from the generalized fluctuation-dissipation relation, and find 
a qualitative agreement between the two different nonequilibrium temperatures. 
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The Fluctuation Theorem (FT) concerns the fluctu- 
ations of the entropy production rate a{t) in the sta- 
tionary nonequilibrium states of a chaotic driven sys- 
tem. We will set /cg = 1 and define cr+ = {cr{t)), where 
(•) is the time average in stationary state; if cr+ > 
the system is out of equilibrium. Defining the variable 
p{t) = {Ta+)~^ J^^"^ ds a{s) (such that {p) — l), its Prob- 
ability Distribution Function (PDF) TTrip), and the large 
deviations function 

CAp) = log7r,(p) , dp) - lim CAp) , (1) 

r — >oo 

the FT states that, for \p\ < p* (where p* is defined by 
limp^ip. ({p) = —00), the following relation -also called 
Fluctuation Relation (FR)- must hold: 

c(p) - a-p) = P<^+ ■ (2) 

The validity of this relation was first shown by Evans et 
al in a numerical simulation of a sheared fluid 1] and 
subsequently proven for reversible Anosov systems by 
Gallavotti and Cohen 2j. Gallavotti then showed that, 
close to equilibrium ((T-|_ — + 0), the FR implies the usual 
Fluctuation-Dissipation Relation (FDR) [Sj. In the re- 
cent past, the FR has been tested under a wide class of 
different conditions, and is now believed to be a very gen- 
eral relation for chaotic systems ^ IE Bi recently, it 
has been also tested in some experiments |3- 

Extending the FT to the case of driven Langevin sys- 
tems, Kurchan pointed out that "the FR might be vi- 
olated for those (infinite) driven systems which in the 
absence of drive have a slow relaxational dynamics that 
does not lead them to equilibrium in finite times" [^. 
This is (by definition) the case of driven glassy systems. 
Driven glassy systems have been widely studied by nu- 
merical simulations: in |^ a uniform velocity gradient 7 
was applied on a Lennard- Jones liquid (which manifests 
glassy behavior below the glass transition temperature 
Tg [iflj) at fixed kinetic temperature T. In presence of 
the driving force, the system becomes stationary also be- 
low Tg, while in the absence of drive the system is not 



able to equilibrate with the bath and ages indefinitely. 
It was shown that below Tg the FDR does not hold any- 
more (because in absence of drive the system is not able 
to reach equilibrium below Tg) but can be generalized, 
for small driving force, introducing an "effective temper- 
ature" Teff -higher than the temperature of the bath- 
associated with the "slow" modes that in absence of drive 
are responsible for the glassy behavior ^| . The break- 
down of the FDR and the close relation between the latter 
and the FR support the conjecture of Kurchan that the 
FR also has to be modified below Tg. 

A possible generalization of the FR, of the form 

Cip) - a-p) - Xpa+ , (3) 

was proposed in Q in the context of chaotic dynamical 
systems. It has also been proposed to define Tpfi = T/X 
as the "temperature" in nonequilibrium steady states 
|l2l | . A similar generalization has been proposed by many 
authors in the context of glassy systems, following the re- 
ported Kurchan's observation, and some attempts have 
been made in order to relate Tpn with the effective tem- 
perature Teff introduced in the generalized FDR p^[l3 |. 
Recently, a connection between the generalized FDR and 
the FR has been derived in a model for the Brownian dif- 
fusion of a particle in a nonequilibrium environment [iSj . 

However, up to now numerical studies of the FR have 
been performed only in the high temperature region 
(T » Tg). The aim of this paper is to test the FR 
below Tg in a numerical simulation of a Lennard- Jones 
like liquid. We measured C,{p) and found that the data 
are consistent with Eq. O with X < 1, while above Tg 
one has X — 1, consistently with what has been found in 
previous works. We measured also the effective temper- 
ature Tf,ff from the generalized FDR and found a good 
agreement between TpR = T/X and T^ff. 

The investigated system is a 80:20 binary mixture of 
N—QQ particles of equal mass m interacting via a Soft 
Sphere Potential (SSP) y„/3(r) = e„/3 (^)^^ {a, (3 G 
[A, B]). This system has been introduced and character- 
ized in equilibrium by De Michele et al as a modifi- 



2 



cation of the standard LJ Kob- Andersen mixture that is 
known to avoid crystaUization on very long time scales, 
and hence to be a very good model of glass former; it has 
been chosen because the SSP can be cut at very short dis- 
tance (~ 1.5a aa) allowing the system to be very small 
(A''=66) in order to observe the negative values oi p that 
are required to test Eq. |21 All the quantities are reported 
in units of m, tAA and oaa- In these units the integra- 
tion step is dt — 0.005. The particles are confined in a 
cubic box with Lees-Edwards boundary conditions [T^ 
at density p = 1.2. The shear flow is applied to the sys- 
tem along the x direction with a gradient velocity field 
along the y axis. The molecular dynamics simulation is 
performed using SLLOD equations of motion [17i |: 



Qi ^ 1- 7qyiX , 

m 



Pi = F^{q) - jpy^x - apr , (4) 



where Fi{q) = ^dq.V{q) and a is a Gaussian ther- 

2 

mostat that fixes the kinetic energy J^- ^ = |iVT. 
The entropy production rate is defined as the dissipated 
power W divided by the kinetic temperature T [l8| : 
a{p, q) = Wip, q)/T = -jP^yip, q)/T, where P.y{p, q) = 
'^i[PxiViii+ qyiFxi{q)] is the xy component of the stress 
tensor [iTj. 
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FIG. 1; Viscosity as a function of temperature for different 
values of 7. The continuous line is a fit to a Vogel-Tamman- 
Fulcher law, r?(T) = rjoexp{:^^) with -qo = 5.2, A = 0.99, 
To = 0.85. ° 



In Fig.^we report the viscosity 77 = {Pxy) H as a func- 
tion of the temperature T for different values of the shear 
rate 7. At 7 = the viscosity seems to diverge at a tem- 
perature To ~ 0.85; however, we are able to equilibrate 
our system only down to T ^ 1.1, that provides an esti- 
mate for the glass transition temperature Tg. For 7 > 
the system becomes stationary and the viscosity is finite 
at all temperatures, even below Tg. 

Very long simulation runs (up to 2 • 10^ time steps) 
have been performed to measure the PDF of the entropy 



production rate along the line 7 = 0.03. During the run, 
p{t) has been measured on subsequent time intervals of 
duration r. From this dataset, we constructed the his- 
tograms of T^rip) and the large deviations function Cr(p) 
defined in Eq. The function C^r (p) is observed to con- 
verge to its asymptotic value ({p) for r > Tq, Tq being the 
relaxation time of the autocorrelation function of crit). 

In the upper panel of Fig. El we report the functions 
Cr(p) for 7 = 0.03 and T = 1.4 > Tg. The asymptotic 
function is obtained for r ^ 5 and can be described 
by a simple Gaussian form, ({p) = —{p — 1)^/2^^, even 
if small non-Gaussian tails are observed. In the lower 
panel of Fig. [21we report Cr(p) — Ct{—p) as a function of 
p(j+. The FR, Eq. 13 predicts the plot to be a straight 
line with slope 1 for large r; this is indeed the case for 
T 5, consistently with what has been found in the 
literature 

In the upper panel of Fig. 13 we report the functions 
Crip) for 7 = 0.03 and T = 0.8 < Tg. In this case, the 
asymptotic regime is reached for t ^ 6; this value is not 
so different from the one obtained in the previous case 
because the change in viscosity (and hence in relaxation 
time) going from T — 1.4 to T = 0.8 is very small at 
this value of 7 (see Fig. ^ . Also in this case the simple 
Gaussian form gives a good description of the data apart 
from the small non-Gaussian tails. In the lower panel 
of Fig. 13 we report Ct(p) — Ct(~p) as a function of pcr+. 
At variance to what happens for T > Tg, in this case 
the asymptotic slope reached for r ^ 6 is smaller than 
1; thus, the FR given by Eq. [3 has to be generalized 
according to Eq. |3 At this temperature, we get X ^ 
0.83 ±0.05. 

In Fig.0I we report the behavior of the violation factor 
X(T, 7 — 0.03) (full circles) as a function of the temper- 
ature T; note that X becomes smaller than unity exactly 
around Tg ~ 1.1, i.e. when the viscosity starts to diverge 
strongly (see Fig. QJ. Below T ^ 0.4, cr+ becomes so 
large that negative fluctuations of p are extremely rare 
and the violation factor is no longer measurable. We can 
conclude that below Tg the FR does not hold, and our 
data are consistent with Eq.|3 where the coefficient X is 
temperature dependent below Tg and equals 1 above Tg. 

Ha ying checked the validity of Eq. |3 following 
Ref. |l2l |. we can define a nonequilibrium temperature 
as Tfr = T/X. Note that if we define a'{t) = W/Tfr = 
Xa{t), the variable p and the function C{p) are not af- 
fected by the rescaling while cr^ = Xa^. Thus, the FR 
for a' is the usual one given by Eq. |3 If one assumes 
that the fluctuations of entropy production rate should 
satisfy Eq.[21 the entropy production rate should be given 
by a'{t) instead of a{t). 

It is interesting to compare the temperature Tpn with 
the effective temperature Tg// that enters the general- 
ized Fluctuation Dissipation Relation. The latter can be 
measured from the relation Teff = D/p, where D is the 
diffusion constant and /i is the mobility of the particles 
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FIG. 2: Top: the large deviation function C,t(p) ~ 
logTTrip) as a function of p for different values of r at 
T — 1.4 > Tg and 7 = 0.03. Error bars are smaller than 
the symbols except on the tails: they are reported only for 
T — 7.5 to avoid confusion. The line is a Gaussian fit to the 
data with r > 5 for p £ [0,2]. Bottom: Ct(p) — Ct-(— p) as a 
function of pcr+. The FR predicts the plot to be a straight 
line with slope 1 (full line) for large r. 



in the considered steady state |9|, ll9j . This relation gen- 
eralizes the usual equilibrium FDR D = /iT; to compute 
the diffusion constant and the mobihty of type-A par- 
ticles we followed the procedure discussed in Ref. [l9l| . 
In Fig. 0] together with X — T/TpR, we report the ra- 
tio T/Tf,ff (open diamonds) as a function of the bath 
temperature T. The two "effective" temperatures have a 
similar qualitative behavior but do not coincide. 

The origin of this discrepancy will be discussed in de- 
tail in ^3 ; roughly speaking, the point is that the modes 
at all frequencies contribute to the entropy production 
rate, while Tg/ / is the temperature of the slowest modes 
in the systems. At the values of 7 we considered, the 
separation between a "fast" and a "slow" relaxation is 
not so sharp: hence, Tpij should be related to an average 
over all the frequencies of the frequency-dependent effec- 
tive temperature. Note also that we are forced to use a 
very small system in order to observe large fluctuations 



FIG. 3: Same plots as in Fig.glfor 7 = 0.03 and T = 0.8 < Tg. 
In the lower panel the dashed line has slope 1 while the full 
line has slope X = 0.83. 

of the entropy production rate, thus, size effects could 
affect the behavior of the investigated quantities. We be- 
lieve however that the qualitative picture is correct even 
if size effects are not completely negligible. Future works 
will hopefully clarify this issue by exploring lower values 
of 7 for which separation of time scales is more marked; 
however, for low values of 7, size effects are more rele- 
vant and the dynamics of the system is very slow; thus, 
very long simulations of bigger systems, requiring a large 
amount of CPU time, are mandatory. 

An interesting microscopical derivation of Eq. |3 was 
proposed by Bonetto and Gallavotti 0|, who related the 
factor X to the dimensionality of the attractive set of the 
system in its phase space. The latter can be measured 
by computing the Lyapunov spectrum, which in this kind 
of system is composed by pairs of conjugated exponents; 
the latter are constructed by pairing the largest exponent 
with the smallest one and so on The prediction of 
01 is that X = P/TV, where V is the number of pairs 
where one exponent is positive and the other is nega- 
tive, and TV is the total number of pairs. If the attractor 
is dense in phase space, V = Af and X = 1. This re- 
lation is very interesting as -if true- it provides a link 
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FIG. 4: The violation factor X = T/Tfr that enters Eq. 
(full circles) and the ratio T/T^ff from the generalized FDR 
(open diamonds) as a function of the bath temperature T for 
7 = 0.03. 

between the effective temperature and properties of the 
phase space of the system. The Lyapunov spectra have 
been measured by the mean of the standard algorithm of 
Benettin et al |2J and are reported in Fig.|Slfor 7 = 0.03, 
T = 1.2 > Tg and T = 0.8 < Tg. Unfortunately, no qual- 
itative change in the spectrum is observed on crossing Tg 
and in particular V/M = 1 above and below Tg. Thus, it 
seems that the theory of does not apply to our model 
below Tg. Note however that this theory is developed un- 
der the assumption of a strong chaoticity of the system, 
while below Tg and for 7 ~ the dynamics of our system 
becomes slower and slower. Thus, our results should not 
be regarded as invalidating the conjecture of 0, but as 
indicating that the hypothesis of [J| (essentially, the re- 
quirement of strong chaoticity) are not satisfied by our 
model below Tg. 

To resume, we studied the fluctuations of entropy pro- 
duction in a numerical simulation of a Lennard- Jones hke 
fluid above and below the glass transition temperature 
Tg. We showed that below Tg the Fluctuation Rela- 
tion does not hold; in particular, our data are consistent 
with a modified form of the PR expressed by Eq. 0| We 
also showed that the behavior of the temperature derived 
from Eq. |31 Tpn, = T/X^ is qualitatively similar to that 
of the effective temperature Tg/ / that is usually defined 
from the generalized Fluctuation Dissipation Relation. 
A relation between TpR and T^ff has been proposed in 
[1^ Q| and our result are consistent with a recent quan- 
titative derivation of this relation in a simplified model 
[l5| . Finally, we tested a conjecture that relates the fac- 
tor X in Eq. O to properties of the phase space of the 
considered system; unfortunately, our data are not con- 
sistent with this conjecture; thus, we believe that the 
violation of the FR is, in our case, of different origin 
than that proposed in \^ . We hope that future work will 



FIG. 5: Lyapunov exponents for 7 = 0.03 and T = 0.8, 1.2. 
For both temperatures each pair consists of one positive and 
one negative exponent, thus V/M = 1. 



clarify this important issue. 
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